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§ : 0. INTRODUCTION. 

o 

■ In our paper we construct a new infinite family of symmetries of the Whitham 

equations (averaged Korteveg-de-Vries equation). In contrast with the ordinary hydro- 
dynamic-type fiows these symmetries are nonhomogeneous (i.e. they act nontrivially at 
the constant solutions), are nonlocal, explicitly depend upon space and time coordinates 
and form a noncommutative algebra, isomorphic to the algebra of the polynomial vector 
^ ■ fields in the complex plane. 

We will study the averaged Korteveg-de-Vries (KdV) equation. But the main results 
^ \ of our paper can be easily extended to other averaged systems, associated with the 
\^ \ integrable by the inverse scattering transform 1-1-1 equations. 
We will assume original KdV to be normalized as: 

13 , , 

nt = -u^xx - -uux- (0.1) 

We will also need a 2 -|- 1 - dimensional generalization of KdV - the Kadomtsev - 
Petviashvili (KP) equation 

1 3 

{ut - -Uxxx + -uux)x = ^a'^Uyy, c? = ±1. (0.2) 

Let us recall the definition of the averaged KdV equations in the classical one-phase 
case (see [1] for more detailes). Equation (0.1) possesses a three-parametric family of 

^This text wes published in the book "Singular Limits of Dispersive Wawes" eds. N.M.Ercolany, 
I.R.Gabitov, C.D.Levermore, D.Serre, Plenum Press, NY, 1994. 
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moving waves 

u{x,t)^2p{x-vt\g2,g3)-v/6. (0.3) 

{p{x, \g2,g3) denotes the Weierstrass p - function). We may try to construct asymp- 
totical KdV solutions of the modulated wave type 

uix, t) = 2pix - viX, T)t + 0(X, T)\g2{X, T),gsiX, T)) - viX, T)/6 + em, (0.4) 

where X = ex, T = et, e is a small parameter (X and T arc called "slow" vari- 
ables). If the correcting term ui is a bounded function asx ~ ^/e,t ~ 1/e then it 
can be proved (see [1], [3]) that the functions g2{X,T), g3(X,T), v{X,T) satisfy some 
hydrodynamical-type equations: 



d 



u , 32iX,T) \ . ( 92{X,T) 

— gsiX,T) \^V{g,,g„v).— \ g,{X,T) (0.5) 



where the matrix V{g2,g3,v) can be expressed via elliptic functions. 

The averaged KdV (0.5) may be transformed to the Riemann diagonal form 



d_ 



/ ri(X,T) \ / v^{f) \ . / n{X,T) 

r2{X,T) \= V2{r) — r2{X,T) |. (0.6) 
\rs{X,T) J \ vs{r) J \ rs{X,T) 



The diagonalizing change of variables (Whitham [1]) is the following : the Riemann 
invariants r^, i = 1,2,3 are the roots of the polynomial R{X) = 4(A -|- — g2{X + 

v/6) — g^, ri <r2 < ra. The matrix V in the diagonal form reads as: 

ri + r2 + r3 2 K{s) 
vi{r} = ^(r2 -ri)- 



3 3' 'K{s)-E{s) 

''^"^ = 3 + 3^^^ - E{s) 

where — (r2 — ri)/(r3 — ri), -E'(.s), K{s) are elliptic integrals. 

Averaging procedure can be treated as a nonlinear analog of the WKB method. 
It can also be applied to the multiphase solutions of soliton equations (nonlinear su- 
perpositions of moving waves) (Flaschka - Forest - McLaughhn [2], Lax - Levermore 
[3]). Averaged KdV equations are known as Whitham equations. Averaged equations 
associated with soliton systems are also known as equations of slow modulations or as 
equations of the soliton lattice dynamics. 

It is well-known that the original KdV equation is an integrable in the Liouville sense 
hamiltonian system. For averaging (at least in the one-phase case) the integrability is 
not necessary but additional structures of the original equations are usually inherited 
in the associated averaged equations. (The connection between the original structures 
and the averaged ones may be very nontrivial). The Whitham equations, associated 
with the KdV solutions with arbitrary number of phases have the following properties. 

1) They can be presented in the hamiltonian form (Dubrovin - Novikov 
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2) They can be written in the Riemann diagonal form (Flaschka - Forest - McLaugh- 
hn [2]). For the hydrodynamical-type systems with more then 2 components the exis- 
tence of the Riemann invariants is a very nontrivial property. 

3) They have an infinite number of conservation laws in involution (i.e. an infinite 
number of mutually commuting symmetries). The full set of integrals consist of 2^1 + 1 
infinite families where g is the number of phases. One of these families is formed 
by the averaged higher KdV flows, the other 2g have no direct analogs in the KdV 
theory. These 2g additional families were constructed by S.P.Tsarev [5] for g = 1 and 
by B.A.Dubrovin [6] for all g. 

4) They can be integrated by the generalized hodograph transform (Tsarev [7], 
[5]). In fact Tsarev proved that any hamiltonian hydrodynamical-type 1-1-1 system 
possessing Riemann invariants is integrable. Algebro-geometrical interpretation of the 
hodograph transform was suggested by I.M.Krichever [8]. 

The fact that the KdV equation has an infinite set of commuting symmetries (higher 
KdV) is well-known. They can be expressed via recursion operator A 

du / A\" d / A\"+^ , 

« = -2x:-T (0-8) 



where A = —d^ + 2d~^udj; + 2u, x = ti, t — t^. But it is less known that the group of 
symmetries of KdV is much wider (the same is valid for all integrable by the inverse 
scattering transform systems). In our parer the averaged KdV symmetries, associated 
with the following KdV symmetries will be studied 

du d f A\"+^ ~ , , / A\\ , , 

as = -'5; (-4) (E((2*+lfe«(-j))'l. (0.9) 

Let all the the higher times tk, k > 3 be equal to zero. Then for n — —2, 0, 2 we 
have: 

du/dp-2^3tu^-2 (Galilean). (0.10) 
1 3 

du/d(3o = ^ti^u^xx - -^uux) + xux + 2u (Scaling). (0.11) 
du/d(32 = ?>t{^u^ ~uu^^^~u^u^^^-^u^u^)^x^^ 

(0.12) 

We call the symmetries (0.9) nonisospectral because they change the spectrum of 
the auxiliary linear problem (see section 2) or Virasoro because they form a noncommu- 
tative algebra, isomorphic to the algebra of the polynomial vector fields in the complex 
plane (the full Virasoro algebra does emerge in the KP theory). The ordinary higher 
KdV are isospectral in this sense. 

First examples of the nonisospectral integrable equations (without connection with 
the symmetry problem) were constructed by F.Calogero - A.Degasperis [9], V.A.Belinskii 
- V.E.Zakharov [10], Maison [11], G.Barucchy - T.Regge [12]. The fact that we can 
generate new KdV symmetries by applying the recursion operator to the scaling trans- 
formation was pointed out by N.H.Ibragimov - A.B.Shabat in [13], but they did not pay 
much attention to this fact because they studied in [13] only local symmetries. Some 
nonisospectral KP symmetries were constructed by H.H.Chen, Y.C.Lee, J.E.Lin [14], 
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F.J.Schwarz [15]. A general approach based on the so-called mastersymmetries was de- 
veloped by B.Fuchssteiner and W.Oevel [16]. A very convenient method of constructing 
nonisospectral symmetries for the integrable by the inverse scattering transform equa- 
tions was developed by A.Yu.Orlov and E.I.Schulman [17]. An analog of formula (0.9) 
for the KP equation was found by A.S.Fokas and P.M.Santini [18]. 

All these papers were dedicated to the algebraic theory of the nonisospectral sym- 
metries or to the vanishing in the infinity case. The periodic (finite-gap) theory needs a 
special consideration. l.M.Krichever and S.P.Novikov calculated the action of a part of 
such symmetries on the finite-gap KP solutions in [19] (their subalgebra was associated 
with a given KP solution and did not vary the spectral curve). They posed also the 
problem for the full KP analog of the algebra (0.9). 

The action of the symmetries (0.9) on the finite-gap KdV solutions was calculated 
by A.Yu.Orlov and the author in [20], [21]. Similar results were also obtained in [20], 
[21] for the KP equation. It was shown that this action has a natural geometrical 
interpretation (see paragraph 1.3). 

Late 80'ies a number of hnks between the nonhnear integerable equations and the 
conformal field theory was discovered. For example the finite-gap r-function in the KP 
theory coincides with the determinant of the d operator on the appropriate bundle (see, 
for example [20]). This determinant plays the central role in the string theory. The 
matrix models of the two-dimensional gravity give us another important example [22] . 
The partition function for the one- matrix model in the double scale limit coincides with 
the r-function, corresponding to a special KdV solution, determined by the constraint 
([22]): 

du 

In the conformal field theory the algebra of the vector fields in the circle and its 
central extension - the Virasoro algebra play the key role. A subalgebra of the algebra of 
nonisospectral symmetries corresponds to them in the KP theory (the central extension 
does emerge if we calculate the action of these symmetries at the r-function) (see [34] , 
[20]) . 

Similar results for the two-dimensional topological quantum theory were obtained 
by l.M.Krichever [23]. He proved that the averaged KP equations are connected with 
the topological models and calculated the action of the Virasoro symmetries on these 
equations in terms of the averaged r-function. Applications of the Whitham equations 
to the topological models were also considered by B.A.Dubrovin [6]. 

In our paper we use the results of [23], [20], [21] to calculate the action of non- 
isospectral symmetries on the averaged KdV in terms of the Riemann invariants. We 
will assume the averaged KdV to be written in the Flaschka - Forest - McLaughlin 
form. The plan of our parer in the following. In the first section we recall the necessary 
definitions and results from the Riemann surfaces theory, including the deformations 
of the complex structures via the vector fields action. In the section 2 we recall the 
necessary results from the periodic KdV and KP theory including the action of the 
nonisospectral symmetries on the finite-gap solutions. The section 3 is dedicated to the 
averaged KdV equations. We recall the construction of the Abel Whitham hierarchy 
and present a new noncommutative set of symmetries. 
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1. THE RIEMANN SURFACES THEORY. SOME DEFINITIONS AND 
RESULTS. 



In this section we recall the constructions from the Riemann surfaces theory we will 
use later ([24] - [26], [6], [20], [33]). 

In our paper the words " Riemann surface" will always denote a compact nondegen- 
erate Riemann surface. Such surfaces may be characterized as close Riemann surfaces 
of finite genus or as algebraic Riemann surfaces. We will consider two main classes - 
the general finite-gap surfaces (they correspond to the KP quasiperiodic solutions) and 
the hyperelliptic Riemann surfaces with a branch point in the infinity (they correspond 
to the quasiperiodic KdV solutions). 

A Riemann surface T is called hyperelliptic if there exists a meromorphic function 
E onT with two simple poles or with one second-order pole. This function maps F to 
the complex plane, the covering : F ^ C is two-sheeted so F is isomorphic to 

a) Y^^{E-E,)...{E-E2g+2)or (l.l.a) 

b) Y^^(E-E^)...{E-E2g+i) (1.1.6) 

respectively (here g is the genus of F). The branch points of F over E are Ei, . . . , E2g+2 
in the case a or . . . , £'2^+1, 00 in the case b. In our paper only the hyperelliptic 
surfaces of the type b will be considered because only they are related to the KdV 
theory. 

1.1. Cycles On The Riemann Surfaces. 

Any compact nondegenerate Riemann surface is topologically equivalent to the 
sphere with a finite number of handles. This number is called genus, we will denote it 
9- 

The canonical basis of 1-cycles in F consists oi 2g elements ai, . . . , a^, 61, . . . , bg 
such, that 

Qi o Qj — bi o bj — 0, Qi o bj — Sij. (1.2) 

where o denotes the intersection number. (Of course there arc infinitely many nonho- 
motopical choices of canonical bases.) The surface F can be described as the result of 
gluing the sides of a 4g(-sided polygon (the picture corresponds to g = 2) 




in the following way: af <-> a~, bf <-> b~ . (All the vertices are glued together). The 

sides af, bf correspond to the cycles Oj, bi respectively. 

Let F be a hyperelliptic surface - a two-sheeted covering over the i?-plane with the 
branch points 00, Ei, E2g+i and the cuts (— oo,£^i), (£'2, £"3), (£25, -£29+1)- 
Then we can choose the following canonical basis 
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where the point £'0 correspond to the vertices of the polygon. 

1.2. Differentials On The Riemann Surfaces. 

Let us have a marked point 00 in our Riemann surface F with a local parameter z. 
We use the notation 00 because this point corresponds to the infinite value of energy in 
the spectral theory. The local parameter 2; is a function defined in the neighbourhood 
of 00 such, that z{(yo) = 0, dz{'~f) |^=oo7^ 0. For convenience we introduce an additional 
function A(7) = l/z{'y). For the hyperelliptic surfaces we will always assume E{oo) — 
00 and E — — A^. All the A-expansions will be defined in the neighbourhood of the 
point 00. 

For us the following objects in F will be necessary. 

1) Canonical basis of holomorphic differentials uui, . . . , ujg. Differentials ujj are 
determined by the normalization conditions 

f uji^ 5ij. (1.3) 

Jaj 

The matrix (Bij) 

Bij = / iJi. (1.4) 
Jbj 

is called Riemann matrix. It is symmetrical (S^ = Bjj), the imaginary part of the 
Riemann matrix Im Bij is positive definite. The coefficients of the Taylor series for cuj 
in 00 we will denote 

= E ^n.A-'-^ciA (1.5) 

m>l 

2) Meromorphic differentials Qj with the only pole in the point 00 such, that 

Qfe = (i(A'=)+0(l), (f ilk^O. (1.6) 

We denote 

/ Qk^{Uk)j, Uk={{Uk)l,...,{Uk)g) 

Jbj 

Qfe = diX") + QkmX-'^-'dX. (1.7) 

m>l 
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The Riemann bilinear relations (see (1.16) below) result in 

{Uj)k = -2mq^j, Qij = Qji. (1.8) 

If r is a hyperelliptic curve then fl2k — d{{—E)''), {U2k)i — 0, Qui = if at least one of 
the indexes A;, / is even. 

3) Multivalued holomorphic differentials a;^, a\ (they are defined only in the hyperel- 
liptic Riemann surfaces). These differentials are determined by the following properties 

0, (1.9) 

0, (1.10) 

(1.11) 

Here A^. and A;,, mean the increment of the differentials when going along the cycles 
Qi and bi respectively. We will denote 

= E Q^A— MA, al=Y: i^LA-^-^A, (1.13) 

m>l m>l 

4) Quasimomentum 1-differential dp, meromorphic in F with the only pole in the 
point oo such, that 

dp — —id{X) + 0(1), Im dp — lm i dp — 0. 

J Qj J bj 

5) Algebra of the holomorphic vector fields in the punctured neighbourhood of 00 
with the basis 

= A^+^9a, -00 < i < 00. (1.14) 
These fields have the following commutators 

[h,lj]^{j-t)li+j. (1.15) 

In the hyperelliptic case we will consider a subalgebra, generated by the elements 
with even indexes. All elements of this subalgebra are single- valued in the £'-plane (in 
the neighbourhood of infinity). 

L, = /2. = 2{-iyE'+'dE. 

The deformations of F, generated by these elements preserve the hyperelliptic struc- 
ture. (The action of the vector fields on the Riemann surfaces will be described below 
in the section 1.3). 

6) Holomorphic 2-differentials in F. These differentials can be written in the local 

coordinates as u^'^^ = uj^'^\z){dzy , the space of such differentials is 3(? — 3 - dimensional 
as (7 > 1, 1 - dimensional as g = 1 and empty as g = 0. We have a standard overdeter- 
mined full set in this space ujIj = uJiUJj where uJi are holomorphic 1-differentials. 

7) Riemann bilinear relations. 



^a.U^l^d.AE"). ^bA 
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Let ^2 be multivalued meromorphic differentials in F such, that 

a) Qi, Q2 have no branch points and are locally single- valued. 

b) All the residues of Qi, Q2 are equal to 0. 

c) d^^Aa.f2j, d~^Afnflj are single-valued meromorphic functions in T without sin- 
gularities on the cycles. We assume all these functions to be equal to in the point of 
intersection of all cycles (this point correspond to the vertices of the polygon) . 

Then we have the following relation: 

2Tii ^ res(d-^l]i)l]2 = ^ { / l^i / - / ^^2 / ^^i + / (d'^Afe, 1^2)^1 1 + 

j w Oj Jbi Jai Jbi J at ) 

+ E|- / {d-'Ab,ni)n2+ (f {d-'AaMn2 - / {d-'AaP2)^i]+ (1.16) 

j I Jai Jbi Jbi ) 

+ lid-'A,,n^){A,,n2) - f^ {d-'Aa,n2){Aa,n,)Y 

All the terms in the right-hand side of (1-16) are correctly defined because of c). 
The property b) guaranties the correctness of the left-hand side. 
Remcirk. All the residues of Qi, Q2 are equal to so we have 

^res(rf-^f7i)n2 = -^res(rf-^02)^^i (1-17) 

8) We will use the following scalar product, introduced in [27] and generalized in 
[23], [6]. Let the 1 - forms fli, O2 be linear combinations of differentials, defined in the 
points 1-3. Then 

Vhif22 = res {d'^Qi)+Q2 - ^ E / ^1 j[ ^2+ 

+ ^E£(^"'^^.^i)^2-^E£(^"'^«.^i)^2. (1.18) 

(We have only one singular point in (1.18) so we omit the J2 sign before the residue). 
The brackets ( )+ in (1.18) denote the singular part of a function 

/m=+oo \ m=+oo 

E = E «-A"^ (1-19) 



+ 



m=l 



Comparing (1.17) with the Riemann relations (1.16) we see that the product (1.18) 
is symmetrical 

Va,n, = Vn,n,. (1.20) 

( The last two terms in (1.16) are equal to 0. ) 
Simple direct calculations result in: 

VhfcQi = Qkl — Qlk 
Vuikiii = Qkl 
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The following property of the scalar product Vh^Qj is important for the averaged 
equations theory. 

Lemma 1.1 Let F be a hyperelliptic Riemann surface with the branch points oo, 
El, E2g+i. Then the scalar products Vh^o,, VL^n,, VL-Hp K«f2,, Vuk'^i, K>fc<^°, VLj^^;" 
can be expressed via branch points Ei, . . . , -E'29+1 and the normalization point Eq in 
terms of the hyperelliptic integrals (elliptic integrals a.s g = 1). 

Proof of the Lemma. The forms 0;^, Qk in the hyperelliptic case have the following 
representations 

{EliAkiE^-')dE 



^-A{E -Ei)...{E- E2g+i) 



^-A{E - El) ... (E - E25+1) 
A; = 2n + l,n > 0. 

The coefficients A^i, Cm are uniquely determined by the normalization conditions 

Ui^Sij, i Qi^O. (1.23) 



Equations (1.23) are equivalent to a linear system on A^i, Cm- All the coefficients of 
this system are hyperelliptic integrals, the system is non-degenerate (see [25]) and we 
can solve it. The coefficients of the Taylor series in 00 for a;^, can be algebraically 
expressed via Aki, Cki, Ej. The remark that all the integrals 

f{E'-E^)uJi, f{E''-E^)ni (1.24) 

are hyperelliptic completes the proof. 

1.3. Deformations Of The Riemann Surfaces And The Riemann Problem. 
The Moduli Space. 
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All the surfaces of genus g are topologically equivalent. But if > 1 then they may 
be different as complex manifolds. Consider the simplest case g — 1 oi complex tori. 
There exists an unique holomorphic differential uji such that 

uji = 1. (1.25) 

■la 

The parameter 

T = j^uji, Im r > (1.26) 

is correctly defined if we have a fixed basis of 1 - cycles a, h in F. The surface F can be 
described as a factor of the complex plane by the group of shifts, generated by: 

z ^ z + l,z ^ z + T. (1.27) 

Let F, F' be complex tori (g=l) with the bases of 1 - cycles a, h and a', h' respectively. 
Then a complex map / : F — > F' such that /(a) = a', f{h) — b' exists if and only if 
t = t'. 

Parameter r depends on the choice of basic cycles. Let us have two bases a, b and 
a, b respectively. Then there exists a matrix ^ ^ ^ e sZ(2, Z) such, that 

f = (1.28) 

The group sl{2, Z) consists of all 2 x 2 integer matrices such that aS — — 1. 
Summing all of this we have: 

Two complex (one-dimensional in the complex sense) tori with the parameters r 
and r' are isomorphic if and only if there exists a matrix 

( " 1 esl(2,Z) such, that r' ^ ^^^J . (1.29) 

The space of all 1-dimensional complex tori is called the moduli space for = 1 or the 
space of all complex structures on the surfaces of genus g = 1- We have proved that 
the moduli space for gr = 1 is the factor of the complex upper-plane Im r > by the 
group sZ(2, Z). 

It is possible to consider the moduli space of all Riemann surfaces with a set of 
marked points and the moduli space of the Riemann surfaces with a set of marked 
points and a local parameter in one of them. The last space is infinite-dimensional. 

The moduli space of the Riemann surfaces with g > 1 is a, 3g — 3-dimensional 
complex manifold. We do not want to discuss its properties in details. But the following 
construction is very important for us. 

Let F be a Riemann surface with a marked point oo, a local parameter z = 1/X 
in the neighbourhood of oo, 5 be a small contour surrounding oo, v = v{X)dX he a 
holomorphic vector field in the vicinity of the contour U{S). Then we can construct 
a family of new Riemann surfaces F^ depending on a parameter a (see [25]). We will 
assume a to be sufficiently small. 

The contour S splits F to a small disk D, containing the point oo and an open 
Riemann surface T\D ( T\D denotes the set of point 7 such, that 7 e F, 7 ^ D). We 
can cover F by two regions F_|_, F_ such, that: 
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1) dg r_. 

2) r\D c r+. 

3) r+ nr_ = u{S). 

If 7 is a point of U {S) then we will denote the correspondent points in r+ and r_ 
7+ and 7_ respectively. F may be treated as the result of gluing r+ to r_ / : r+ — > r_ 
/(7+) = 7_. Let us introduce a new gluing function by formula 

/a(7+) = exp(a't;)7_, (1.30) 

where exp(Q;t;)7_ denotes the shift of the point 7_ via the vector field v after the lapse 
of the time a. Then we can define a new Riemann surface Fq, as the result of gluing 
r_|_ to r_ via the function f^. 

The Riemann surfaces F and are constructed from the same parts F_|_ and F_. 
Thus we have a map E : F — > Fq, coinciding with identical maps F_,_ — > F_,_ on r\D 
and F_ F_ on D. Map E has a jump on the contour S. If we have some marked 
points in F or local parameters in some of them then we will map these objects by E. 

Let us have a infinitesimal transformation (a is infinitely small) and A be a holo- 
morphic tensor field in F^. The map E carries A from F^ to F so we can treat A as a 
tensor field in F with a jump on the contour S. The boundary values of A on 5" A+ 
and A_ satisfy the following relation 

A+-A_ = aL^A, (1.31) 

where L„ denotes the Lie derivative 

Lvwa, 9{\){dXr = {v{X)g\\) + av'{X)g{X)) (rfA)" (1.32) 

(a is infinitely small so we can write A_|_ in the right-hand side of (1.31) as well as A_). 

Formula (1.31) gives us a very convenient method for calculation the variations of 
tensor object via complex structure variations: we calculate the right-hand side of (1.31) 
and solve the Riemann problem in the appropriate functional class. Riemann problem 
is one of the basic objects in the soliton theory so this approach allows us to connect 
these deformations of the Riemann surfaces with the non-isospectral symmetries (see 
section 2.2). 

In the hyperelliptic case our construction has the following interpretation. Let F 
be a hyperelliptic surface over the plane with the branch points oo, Ei, . . . , E2g+i, 
V — v{E)dE be a vector field in the £'-Plane, nonsingular for all E ^ oo, the local 
parameter in the point oobe^ = l/A, A = \J —E. 

Then Fq, is a hyperelliptic Riemann surface with the branch points oo, exp(af)i?i, 
exp(Q;t')ii^2s+i and the local parameter z = 1/A, A = ^ —E. The map E carries 
the point E± in F to the point exp(Q;i;)ii^± in F^ if the point E± is located outside the 
neighbourhood of oo surrounded by a small contour S and carries the point E± in F to 
the point E± in Fq if the point E± is located inside the neighbourhood of oo (the sign 
± means the upper or lower sheet). 

If we have a family of deformed hyperelliptic Riemann surfaces Fq, and a function 
on this family f{E±, a) then we have two differential operators 

9 J = 1^ (1-33) 



11 



and 

df 

DJ=^ + LJ, (1.34) 

where is the Lie derivative (1-32) in the first argument. In contrast with da the 

operator can be generahzed to arbitrary Riemann surfaces and can be treated as a 

connection generated by the map E. 

The following property of the connection Da is very important for us: 

Lemma 1.2. Let A{E±,a) be a holomorphic tensor field in rQ,\oo for all a. Then 

DaA{E±,a) is nonsingular outside the neighbourhood of oo (for da in is not valid in 

the branch points). 

While calculating the nonisospectral symmetries we will need the derivatives of the 
holomorphic differentials by the complex structures. In the hyperelliptic case wc have: 

Lemma 1.3. Let F be a hyperelliptic Riemann surface, L„ = = 2(— 1)"£'"+^(9£; 
be a holomorphic vector field in the E-pla.nc, i > — 1, Fq be the deformation of F via 
the vector field v. Then 

2n-l 

Da^k — kD,2n+k + X/ Qkm^2n-m, (1.35) 
m=l 

2n-l 

Dal^k = J2 (lkm^2n-m: (1-36) 
m=l 

DaUjt = 2(-l)"/c<+, + QlAn-m. (1-37) 

m=l 

DaOt = 2(-l)-A;<+, + RtJl2n-m. (1-38) 

m=l 

Proof of the Lemma 1.3. In accord with the rule described above we have to solve 
the following Riemann problem: 

Let A be one of the differentials fife, a;^, o;^, o"^. Then we have to construct a pair 
of differentials A+, A_ such that 

1) A+ - A_ = L^A. 

2) A_ is defined and nonsingular in the neighbourhood of oo. 

3) A_|_ is nonsingular for all E± ^ oo. 

4) 4 ^+ = 0- 

5) A+ is single- valued as A = Qfe or A = Uk, 



AaA+ = 2(-l)"A;A,^<+fe, A,,A+ = as Ah 



UJi. 



A,A+ = 2(-l)-A;A„,<+„ A„,A+ = as A+ = a," 
Then DaA = A+. 

In fact wc arc looking for a differential with the properties 3-5 in the finite part 
of F and the prescribed singularity L^A in oo. But such differential can be easily 
constructed as a linear combination of flk, Wk, '^ti '^t- This completes the proof. 

Corollary 1. ([25]). 

^ (1.39) 
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where the product ^(A)^^ uJi{X)dX ujj{X)dX is the 1-form v{X)uJi{X)ujj{X)dX. This for- 
mula is vahd for general Riemann surfaces as well as for hyperelliptic ones. We assume 
that the contour S goes around the point oo counterclockwise. 

Corollary 2. The vector field v do not vary the complex structure of T if and only 
if all the integrals 




is are equal to 0. Here ool is the basis of holomorphic 2-forms. It proves that the 
moduli space is 3gi — 3 dimensional as 51 > 1 and 1 - dimensional as 51 = 1. 

1.4. Riemann Theta- Functions. 

Let bij be a complex symmetrical g x g matrix such, that Re bij is negative definite, 
zhe a complex g^-component vector. Then the Riemann theta-function can be defined 
as an infinite sum ([26]) 




(1.40) 



where m^, k — 1, . . . , g are arbitrary integers. This sum converges for all z. The theta 
function has the following periodicity properties 

9{zi,Z2, ...,Zk + 27ii, . . .,Zg\bij) = 9{zi,Z2, ...,Zk,.. ■,Zg\bij). (1.41) 

d{zi + bik, Z2 + b2k, ■■■,Zg + bgk\bij) = 9{zi, 2^2, ... , Zg\bij) exp{-&fejfc/2 - Zk}. (1-42) 

The zeros ot the theta function are described by the following Lemma (see [26]) 
Lemma 1.4. Let F be a Riemann surface with a marked point 70, Bij be the matrix 

of periods (1.4), bji = 2-KiBji, A(j) be the Abel transform, i.e. ^(7) is a multivalued 

map r — > determined by the formula 

{A{^))k = 27ri r ujk, k^l,...,g. (1.43) 

Then there exists a g'-dimensional complex vector K — X(r,7o) such that the 
function 

(^(7) = 9{A{^) - l(7i) - ... - A(7,) + K\kj), 7 e r. (1.44) 

has one of the following properties: 

1) (^9(7) = or 

2) (^(7) has exactly g zeros in the points 71, . . . , ■jg. 

— * 

K is called Riemann constants vector. 

2. PERIODIC THEORY OF THE KORTEVEG - DE - VRIES AND 
KADOMTSEV - PETVIASHVILI EQUATIONS. 

In this section we will recall some facts from the KdV and KP theory. Two topics 
are the most interesting for us: periodic (quasiperiodic) finite-gap theory and the action 
of the nonisospectral symmetries on the finite-gap solutions. 
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Finite-gap KdV solutions can be treated as nonlinear superpositions of the moving 
waves. Such solutions are the basic objects for the averaging procedure. The have 
been constructed in the papers of S.P.Novikov, B.A.Dubrovin, V.B.Matveev, A.R.Its, 
P.Lax, H.McKean and P. van Moerbeke (see book [28] for more detailed description and 
references). Finite-gap KP solutions were constructed by I. M.Krichever [29]. The direct 
periodic problem for KP is more complicated (the results, obtained by I. M.Krichever 
can be found in [29]). 

The first step in constructing averaged nonisospectral symmetries is the following: 
the action of these symmetries on the finite-gap solution is calculated. This problem was 
solved by A.Yu.Orlov and the author in [20], [21]. We recall some results of these parers. 
We use the representation for nonisospectral symmetries, suggested by A.Yu.Orlov and 
E.I.Schulman [17] as the most convenient for us. 

2.1. KdV And KP Theory. Integration And Isospectral Symmetries. Peri- 
odic Theory. Bciker-Akhiezer Function And Cauchy Kernel. 

The theory of the KdV equation 

ut = -^u^xx - ^uu^- (2.1) 
is based upon the existence of the following representation. Let 

L^-dl + u{x, t), A^d^^- ^{ud^ + d^u) (2.2) 

be ordinary differential operators depending on an extra parameter t. Then the function 
u{x,T) satisfy (2.1) if and only if the following relation takes place 

dL/dt = [L,A]. (2.3) 

Representation (2.3) is called Lax pair or L — A pair for KdV. 

One of the first results of the soliton theory was the existence of infinitely many 
mutually commuting KdV symmetries. They can be written via the recursion operator 

u = -2—(--] -1, (2.4) 



dt2n+\ dx\Aj dx\ A 
where 

A = -9^ + 2d-^ud^ + 2u (2.5) 

or in the Lax form 

9V9w = [^>^2n+l], (2.6) 

where 

^2n-M={(-i^)'^}+, (2.7) 

(— L)^^ denotes a formal pseudodifferential operator, i.e a series in dx with finite 
number of positive terms and infinite number of negative, { }+ denotes the differential 
part (see [31]). We mark the times by odd indexes to have unified notations for KdV 
and KP. Here ti — x, — t. 
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Equations (2.3), (2.6) result in the following property: the spectrum of L does not 

depend upon the times t^, t^, t2k+i, This is the reason why we call these 

symmetries isospectral. 

KdV equation is integrated by the inverse scattering transform, i.e. we consider the 
"scattering data" for L as a new variable instead of u{x). (We write "scattering data" 
in quotation marks to stress that these data coincides with the physical scattering data 
only for some functional classes of potentials). Because of the isospectral property the 
evolution law for the "scattering data" is very simple (sec [28]). The map from u to 
the scattering data for the vanishing in the infinity potentials in the small amplitude 
limit coincides with the Fourier transform and can be treated as its nonlinear analog. 

The realization of this scheme depends on the functional class of the potential. Let 
us recall the scheme for the periodic case 

u{x + U,t) ^u{x,t), (2.8) 

where u{x,t) is a real nonsingular potential. 

We consider the following spectral problems for L: 

a) Main problem Lil) — Eil^, if) is bounded in x. 

b) Auxiliary problem Lip a = Eipa, V'a(O) = 'ipaO^) = 0. 

The spectrum of the main problem consists of a set of intervals [Ei,E2], [E^jE^], 
[E2n-i, E2n], ■ ■ ■ where Ei < E2 < E3 < E4^ < E^ . . ., E2n+i - £^2n ^ as n ^ 00. 

The spectrum of the auxiliary problem consists of an infinite number of points 
di < d2 < ds < . . . located in the gaps di e [E2, E3], ^2 e [E4, E^], ds e [Eq, E7], . . .. 

The Bloch eigenfunction ip{x. E) normalized by the conditions iIj{x + Il,E) = 
exp{Ilip{E)) iIj{x,E) and iIj{0,E) = 1 is meromorphic on a two-sheeted Riemann sur- 
face r over the E-plane with branch points Ei, E2,. . . ,00, and has simple poles in the 
points 7i, 72, . . ., such that the projection of 7„ to the E-plane coincides with d„. 

The function p{E±) is defined in F and is called quasimomentum. 

The spectrum corresponding to a general potential has infinite number of gaps. But 
for us the so-called finite-gap case when i?2n = -E'2n+i for all n > g is the most impor- 
tant. (The infinite-genus case can also be studied [32] but the answers are much more 
complicated). Finite-gap potentials are dence in the space of all periodic potentials. 

The inverse problem data in the finite-gap case is the following: 

1) 2g + 1 real numbers Ei, E2, . . . , -£'29+1, Ei < E2 < . . . < -£'25+1. The points E^ 
are the boundary point of the spectrum. 

2) g points 71, . . . , 7^ in a hyperelliptic Riemann surface F with the branch points 
00, El, E2, . . . , E2g+i, such, that -£(7^) G [E2k, -£'2fc+i] where E is the projection to the 
£^-plane. 

This data uniquely determines the potential u{x). 

The flows (2.3), (2.6) do not move the branch points E^ but they shift the divisor 
7fc. The evolution of the points 7^ can be described in terms of ordinary differential 
equations, derived by B.A.Dubrovin (see [28]). This system is nonlinear, but it has a 
very nice explicit solution. 

— * 

Lemma 2.1. (see [28]). Let ^(7) be the Abel transform defined in the paragraph 
1.4. Then 

l(7i) + 1(72) + . . . + A{-fg) = Ao + xUi + m + hU5 + ..., (2.9) 
where is some constant vector, the vectors Uk are defined by (1.7). 
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Inverse transform to (2.9) can be written in terms of the thcta- functions. The 
answer is given by the A.R.Its - V.B.Matveev formula (see [28], [33]). 

u{x, h, h, ...) = -2dl log^(Vo(7i, ...,lg) + xU, + tUs + UA, + --- + C(r), (2.10) 

where Vo(7i, . . . , 7g), C(r) are some constants, — 2TriBij, B^j is the Riemann matrix 
(1.4). 

General finite-gap solutions are quasiperiodic. The characterization of periodic so- 
lutions in terms of the inverse data is a complicated problem. 

A slightly different approach to the inverse problem is more convenient in some 
situations (see [28]). Instead of normalizing ip{E, 0) = 1 for all t we consider a function 
^(7, t), 7 G r, t = (ti, t3, ts, . . .), ti = X, ts = t such that 

1) L*(7,i)-^(7)vl>(7,t). 

2) *(7, X + n, ^3, ^5, . . .) = exp(inp(7))*(7, x, ^3, ig, ■ ■ ■) 

3) ^(7,0,0,0,...) = 1 

4) (<9t2„+i — ^2ra+i)^(7,^) = 0. (Symmetries (2.6) is mutually commuting so the 
condition 4 is self-consistent). 

^(7,t) is called Baker- Akhiezer function. It has the following analytical properties: 
Pr.l) ^(7, is meromorphic in r\oo. 

Pr.2) ^(7, has simple poles in the points 71, . . . , 7g and no other singularities in 

r\oo. 

Pr.3) ^^(7,^) has an essential singularity as 7 — > 00 

*(A, i) = exp(e(7, i))[l + Xi{i)/X + X2{i)/X' + . . .] (2.11) 

where in the correspondent to KdV hyperelliptic case A = ©(7, t) = Ax + XH + 

X% + .... 

Lemma 2.2. The properties Pr.l -Pr.3 uniquely determined the function ^^(7,^). 
It can be expressed in terms of the theta- functions (see, for example review [33]) 

^ . ^ e{j:hu, + 1(7) - i(7i) - • ■ ■ - ^(7,) + K)ei-A{^,) - ... - a{^,) + k) 

OiEtkUk - l(7i) - ■ ■ ■ - ^(7.) + K)e{A{^) - l(7i) - ... - Ailg) + K) 

•exp(X]ife f^lk). (2.12) 

— * 

Here K is the vector of Riemann constants (see paragraph 1.4), the integrals are nor- 
malized by P Q,k^ X^ -\- 0(1). 
Comparing (2.12) with 

u{r)^2d,Xi{^ (2.13) 

wc can easily derive (2.12). 

Now we will recall the finite-gap inverse scattering transform for KP [29]. 

1 3 

{Ut - -Uxxx + -^UUx)x = ^Uyy. (2.14) 

Let r be arbitrary Riemann surface of genus g with a marked point 00, a local parameter 
1/A in 00 and g marked points 71, . . . , 7^. Then for the data of general position there 
exists a unique function ^(7,t), 7 G F, t = (ti, ^2, ^3, ^5, • • •), ^1 = x, t2 = y, = t 
with the analytic properties Pr.l - Pr.3 (for KP 6(7, t) = Xx + X^y + XH + XH4 + X% + 
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. . .). Representation (2.12) is valid for general Riemann surfaces but we have to sum 

in (2.12) by even indexes as well as by odd ones. 

For all k there exists a unique ordinary differential operator such that 

(9i„ - i„)*(7,i) = 0, A^^d: + ..., A,^-L. (2.15) 

The compatibility conditions 

[(dt,-A,),{dt,-A,)]^0 (2.16) 

give the Lax pair for KP as A; = 2 and isospectral symmetries as A; > 3. The potential 
u(i) is given by shghtly changed (2.10) - the sum is over all indexes - odd and even, y 
- independent KP solutions satisfy KdV. 

We need also the Baker- Akhiezer conjugate differential ^"'"(7,^). Its analytical 
properties are the following 

Pr.l') ^+(7,t) is a holomorphic in r\oo 1-differential. 

Pr.2') ^"*"(7, has simple zeros in the points 71, . . . , 7^. 

Pr.3') ^+(7, has an essential singularity as 7 — > 00 

*(A, t) = dX exp(-0(7, t))[l + (t)/A + xtii)/)^' + • • •] (2-17) 

where 6(7, i) ^ Xx + X^y + XH + XH4 + X% + . . .. 

(9,„+i:)*+(7,t) = 0, (2.18) 

where ^4+ is the formal conjugate to An, i.e. {a{x) (dx)"')'^ = {—dx)"'a{x). The following 
ortogonality properties are important: 

Lemma 2.3. 1) Let 5" be a small contour surrounding the point 00. Then 

*(7,i>+(7,i')=0 (2.19) 
for all t, F. 

2) Let p{'j) = p dp where dp is the quasimomentum differential defined in the 
section 1, ^(7) be a contour in F consisting of all points 7' such that Im ^(7') = Im ^(7) . 
Then in the contour G(7) the following relation is vahd 

/' + OC 
^(7, t)^+(7', t) = 27ri5(7 - 7'). (2.20) 
-00 

We have to calculate deformations of the Baker - Akhiezer function. For this purpose 
we need an appropriate Cauchy kernel. 

Lemma 2.4. ([20]). The Cauchy - Baker - Akhiezer kernel a;(7,7', with the 
following properties: 

1) ci;(7, 7', t) is a function in 7 and a 1-form in 7'. 

2) a;(7, 7', t) is meromorphic function of 7 in F\oo with simple poles 71, . . . , 7^,, 7'. 

3) As a function of 7' a; (7, 7', is meromorphic in F\oo with one pole 7 and zeros 
in the points 71, . . . , 73. 

4) c<j(7,7',t) = o(exp(9(7, t))) as 7 ^ 00. 
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can be written in the following form 

a;(7, y, x,y,t,...)^^ f ^(7, x', y,t,.. .)*+(7', x', y,t,.. .)dx' . (2.21) 

If 7' ^ ^^(7) the sing in limit of integration in (2.21) is uniquely determined by the 
convergence condition. For 7' e ^(7) the property (2.20) guaranties correctness. 

First formula similar to (2.21) was obtained by I.M.Krichever and S.P.Novikov in 
[19] for systems with discrete x. 

2.2. Nonisospectral Symmetries. Action On The Finite-Gap Solutions. 

We have pointed out in the introduction that KdV equation possesses the following 
set of symmetries 

a;- = -2^(-l) (E((2* + l)t«.(-j))-l. (2.22) 

They can be written in much more convenient form suggested by A.Yu.Orlov and 
E.I.Schulman [17]. 

-2# res U=oo (A^-+^9,*(A,t))*+(A,t), (2.23) 



dr^n dx 

Ordinary higher KdV equations (2.4) have similar representation 

-2^- res U=oo A^"+^^(A, t>+(A, t). (2.24) 



The formulas (2.22), (2.24) can be treated in the following way. Let us substitute the 
asymptotical expansions (2.11), (2.17) to the linear problem 

L*(A, = -A^*(A, i), L*+(A, i) = -A^*+(A, i). (2.25) 

From (2.25) we can calculate all the coefficients Xk{t): xti't) via u{i) (in nonlocal form). 
Then we substitute them to (2.23), (2.25). All the exponents in are reduced and 
the residue can be explicitly calculated. We obtain a close system on ii(i) (may be 
nonlocal). 

Lemma 2.5. Equations (2.23), (2.25) coincides with (2.22), (2.4) respectively. 
Proof of the Lemma. For small n we can check it by direct calculations. Then we 
apply the following identity 

A(*(A, t)*+(A, i)) = -4A2(^(A, t)*+(A, i)). (2.26) 

(it is a direct consequence of (2.25)). 

For the KP equation we can write two-parametric set of symmetries ([17]) 

=-2-^resU=oo(A"ar^(A,t))^+(A,t) (2.27) 



^^nm dx 

but only the symmetries with m — 0,1 are compatible with the finite-gap structure. 
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Theorem 2.1. Let F be a Riemann surface with a marked point oo, a local param- 
eter 1/A in oo, a set of points 71, . . . , 7^ in F where g is the genus of F, /„ = A^+^^a be a 
vector field in the punctured neighbourhood of 00. Consider the deformation of F gen- 
erated by the field (it was described in the paragraph 1.3). Then the correspondent 
variation of the KP solution, constructed by this data reads as 

^ = -2^ res U=oo (A"+^9a*(A, i))*+(A, i) (2.28) 

i.e. it coincides with a symmetry (2.27) such that m = 1. 

Proof of the theorem. In the paragraph 1.3 we have explained that the calculation of 
the Baker-Akhiezer function variation is equivalent to the following Riemann problem 
on the contour surrounding the point 00 

(5^(A,t))+ - (5^(A,t))_ = A"+^aA*(A,t). (2.29) 

Solution of (2.29) reads as 

S-^iX.t) = <f u;{X,i^,i)i^^+'d^^{i^,e). (2.30) 

s 

where uj{X,fi,i) is the Cauchy - Baker - Akhiezer kernel defined in the Lemma 2.4. 
Expanding (2.30) as A — > 00 and using (2.20) we obtain (2.28). 

Corollary 1. Symmetries (2.23) act on the finite-gap KdV solutions as 

dEjdr^n = 2{-irE:+\ dE{^,)/dT2n = 2{-irE^+\^k). (2.31) 

i.e. all the spectral data is shifted via the vector field 2(— 1)'*£'"+^9/9£' (see paragraph 
1.3). 



3. WHITHAM EQUATIONS. THE FULL ABEL HIERARCHY AND 
NONISOSPECTRAL SYMMETRIES. 

The KdV equation (2.3) and the symmetries (2.6) (they are called higher KdV) 
form a commutative set of flows. This set is called KdV hierarchy. 

Averaged KdV hierarchy was constructed in [2], [3]. We will not discuss how the 
averaged KdV equations can be derived from the original ones and so we will only recall 
the answer. 

The starting point for the averaging procedure is the space of all 5^-gap KdV solu- 
tions. Such solutions are parameterized by the branch points E-i, . . . , £^2^+1 and the 
points 7i, . . . , 7g in F. The branch points Ek are integrals of motion and the points 7^ 
play the role of phases (see for example [28] and references therein). 

If we consider a slow modulated wave-type solution then the points E^ slowly depend 
on coordinate and times E^ = Ek{X, T, T5, . . .) where X — ex, T — et, = et^, . . . , 
e ^ 1. Functions X, T, T5, . . . are called slow variables. 

For averaging any full set of integrals can be used. But direct calculations for 
g = 1 (see [1]) shows that the variables Ek result in the simplest form of the averaged 
equations so it is very natural to use them for higher genera. 
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The averaged KdV hierarchy can be written in the following form, suggested by 
Flaschka - Forest - McLaughhn [2] 

Wk [hi,. . . ,h2g+l)—-, [6.1) 



dT2n+l ' V^.,---,^^9+Vq^, 

where 

Wfc ■ ■ ■ ,£/2g+lj - Q^^^^^ ' l"^-^) 

flk are the differentials defined in the paragraph 1.2. We see that the flows (3.1) have 
the Riemann diagonal form. 

The averaged KdV equations have wider symmetry group then the original KdV. 
For example the scaling transform X — > aX, T aT has no analogs in original 
equations. 

From the results of S.P.Tsarev [5], [7] it was known that the Whitham equations 
have 2g + 1 infinite series of symmetries but the averaged KdV hierarchy gives us only 
one of them. The full set of symmetries was constructed by S.P.Tsarev in [5] for g = I 
and B.A.Dubrovin for all g. It contains 2g +1 infinite scries plus one finite and can be 
written as 

dT-^ ^i{Eu)dX' 
dEk _ ajXEk) dEk 

art ni{Ek) ox' ^ ' ' 

dEk _ UjjEk) dEk 
dTl' ~ Q,{Ek) dX' 

where i — 1, . . . , g, n — 1,2, Let us denote 

(see the definitions in the paragraph 1.2). 

We have four families of differentials, times and velocities. To avoid too long nota- 
tions we will use the following agreement: 

may denote any of the differentials Qk, ^i, <^k'i "^f " '^a ^^c corre- 
spondent velocities and times, = ^ngn^ where V is the scalar product defined in 
the paragraph 1.2. 

The following statements are important for us. 

Lemma 3.1 ([2], [6]). The flows (3.1), (3.3) can be written as 



Here we always assume 



dT^ ~ dX ■ ^"^-^^ 



d _ _d_ 
gfF-gfF lE=const 



Proof of the Lemma. The differentials in the both sides of (3.5) have the following 
properties. 
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1) They are single- valued in F. 

2) Their integrals by flj-cycles are equal to 0. 

3) They have second-order poles in the branch points Ei, . . . , £^2g+i said no other 
singularities. All their residues are equal to zero. 

4) Let z — y/E — Ek be local parameter in the neighbourhood of the branch point 
Ek. Then the singular part of both differentials is equal to Q^{z)/2z'^. 

Properties 1-4 uniquely determined a differential so the left-hand side is equal to 
the right-hand one. 

Lemma 3.2 ([2], [6]). All the symmetries (3.1), (3.3) are mutually commuting. 

Proof of the Lemma. Consider the differentials dVt^/dT^ , dVL^ /dT^ . They satisfy 
the properties 1-3 of the Lemma 3.1 and have the same singularities in the branch point 
so they are equal 



and we have 



_^ d_ 

dTfdTj 



d dnf, d dn^^ d d 



dXdTE 



dXdT^ 



Lemma 3.3 ([8], [23], [6]). The scalar products Qfg satisfy the following relation 



07 



(3.6) 



Thus we can define functions 



(3.7) 



Here the vector contains all the times. 

Proof of the Lemma. Prom the formulas (1.21) we see that the coefficients dQ^^/ dT^ 
can be expressed as integrals or expansion coefficients for the differential dfl^ / dT^ Ap- 
plying the Lemma 3.2 we complete the proof. 

Lemma 3.4 Let us shift the branch points Eg of the surface P and the normalization 
point £^0 via a vector field l2k = '2{—^)^Ek+idE 



dE., 



dr. 



2n 



2(-i)"£;: 



n 77in+l 



s = 0,1,..., 2^ + 1. 



(3.8) 



Then we have 



a„,,2k+l n-l 

^^^^ = (2A;+l)w;f+2"+i-wf 

C''^2n rn=0 



Hk 



dT2n 



n-l 

+ E 

m=0 




— W 



2k+l 



Mo 



W 



2n-2m-l 



= 2(-l)'*A;<'=+" 



^2n+l^aik 



n-l 

+ E 

m=0 



d 



a,k_^_\ 

' dX 



Q 



2m+l 



W 



2n-2m-l 

(3.9) 
(3.10) 



w 



2n-2m-l 



(3.11) 
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^TT^ = 2(-l)"fc«;^'=+" 

dT2n 



n-l 

■E 

m=0 L 



d 



d 



.on 



w 



,2n-2m-l 



(3.12) 

(The functions Qj were defined in Lemma 3.3). 

The Lemma is proved by direct calculations using (1.35) - (1.38) and (3.2), (3.4). 
Lemma 3.5. Let us have a pair of flows 



dr 



Rs{Es)+(T.UT)wf^^ 



dE, 



dE, 



w 



FydE^ 

dX' 



(3.13) 



(3.14) 



Then the flows (3.13) and (3.14) commute if and only if the following compatibility 
condition holds: 



dr 



d 



F-y 



-] 

dX 



w 



Fa 



where dw,"' / dr denotes the variations of velocities via the shift 



dE^ 
dr 



Rs{Es). 



This Lemma is proved by simple direct calculation. 

Now wc have prepared everything to formulate and prove our main result. 
Theorem 3.1. The flows 



(3.15) 



(3.16) 



1^ = 2(-1)"E;+i 

dT2n 



Y,{2k + l)T2k+iW 



2k+2n+l I dE, 



. k>0 



dX 



+ <'2(-l)"$:A;(W''^" + rNf"'")+ E Q^k+iw 



k>0 



0<k<n 



2n-2fc-l dE, 

■ dX' 



(3.17) 



commute with the whole Whitham hierarchy (3.1), (3.3). The nonlocal functions Qj 
were defined in the Lemma 3.3. All the partial derivatives of Qj can be expressed via 
branch points E, and normalization point Eq in terms of hyperelliptic integrals (see 
Lemma 1.1). 

Comparing the Lemmas 3.4 and 3.5 we prove this theorem. 



The author is very grateful to B.A.Dubrovin, I.M.Krichever and S.P.Tsarev for 
useful discussions. 
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